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Abstract. Convectively driven dynamos with rotation generating magnetic fields on scales large 
compared with the scale of the turbulent eddies are being reviewed. It is argued that such fields 
can be understood as the result of an a effect. Simulations in Cartesian domains show that such 
large-scale magnetic fields saturate on a time scale compatible with the resistive one, suggesting 
that the magnitude of the a effect is here still constrained by approximate magnetic helicity 
conservation. It is argued that, in the absence of shear and/or any other known large-scale 
dynamo effects, these simulations prove the existence of turbulent a 2 -type dynamos. Finally, 
recent results are discussed in the context of solar and stellar dynamos. 
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1. Introduction 

Stars with outer convection zones are known to display magnetic activity, often in a 
cyclic fashion like in the Sun. Such activity can generally be explained by a turbulent 
dynamo influenced by rotation and stratification to produce the anisotropies required 
for the gene ration of large-scale magn e tic fields. The basic th eory is reasonably well 
understood (|Moffatdll978t |Parkerlll979t iKrause fc Radle"nll98fl . but there continues to 
be substantial controversy until the present day. A major stu mbling block has been 
the understanding of what is kn own as catastrophic quenching (IVainshtein fe Cattaneo 
1992 ; Cattaneo fc Hughes 1996 ) and resistively limited saturation ( Brandenburg! 200 ll) . 
as well as the very existence of the a effect in con vection even without nonlinearity 
( Cattaneo fc Hughed[200l iHughes fc Cattaneoll2008l) . 

The first two issues have been reviewed in detail by IBrandenburg fc Subramanianl 
(2005). The purpose of the present paper is to review recent progress concerning con- 
vective dynamos. However, in view of applications to solar and stellar dynamos, it is 
important to realize that we are still lacking simulations that reproduce the salient fea- 
tures of the solar dynamo. We should therefore keep our eyes open for new phenomena 
that may emerge as simulations become more realistic. 



2. Excitation conditions of small-scale and large-scale dynamos 

Small-scale and large-scale dynamos are quite different in nature. The difference be- 
comes most evident in the nonlinearly saturated regime, provided one allows for what 
we call scale separation, which means that the size of the domain is large compared with 
the scale of the largest (energy-carrying) eddies of the turbulence. In Figure Q] we show 
spectra highlighting the remarkable difference between the two types of dynamos. Con- 
versely, if there is insufficient scale separation, a large-scale dynamo becomes impossible 
and both typ e s of s imulations would look very similar, as has been demonstrated by 
lHaugen et all (|2004 : see their Fig. 23. 

There can be different types of large-scale dynamos. The most frequently studied ones 
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Figure 1. Kinetic and magnetic energy spectra in a turbulence simulation without net helicity 
(left) and with net helicity (right) for a magnetic Prandtl number of unity and a mesh size 
is 512 3 meshp oints. Notice the pronounced peak o f M( k) at fc = fci i n the case with helicity. 
Adapted from Brandenburg fc Subrama nian (2005|) and [Brandenburg (2009), respectively. 

are the a 2 and afl type dynamos. These are dynamos that can produce large-scale 
magnetic fields owing to the presence of kinetic helicity in the turbulence, giving rise to 
an a effect. The presence of shear can further modify the dynamo, making it usually 
easier to excite and favoring oscillatory over non-oscillatory solutions. 

Shear can be a typ i cal re s ult of rotation of a gaseous body in the presence of anisotropic 
turbulence (|Rudigerlll980L Il989f ) . Shear alon e is often found to give rise to large-sca le 
fields - even if the turbu lence is non- helical ( Brandenburg! 2005 ; Yousef et al. 2008al lb: 
Brandenburg et al. 20081) . The nature of such dynamo action is still a matter of debate 



and ranges from incoheren t afl dynamos Jvishniac fc Brandenburg! 1997 : Proctor 2007) 
to shear-current dynamos ( Rogachevskii fc Kleeorinf 20031 l2004j). 

Returning to the a 2 and afl dynamos, it is important to realize that their excitation 
conditions are generally quite different. The onset of small-scale dynamo action depends 
generally on the value of the magnetic Reynolds number, 



Re 



M 



(2.1) 



where M rm s is the typical rms velocity of the turbulence, rj is the microscopic magnetic 
diffusivity, and kf is the forcing or integral wavenumber, i.e. the wavenumber of the 
energy-carrying motions. This is roughly where the peak of the energy spectrum is lo- 
cated. The critical value, -R m ,crit, above which dynamo action commences, depends on 
the value of the magnetic Prandtl number, Pt m = v/v, where y is the microscopic kine- 
matic viscosity and is about 35 for PrM = 1 (jNovikov et al. 1983 : Subramanian 19991 : 
Haugen et al.l 120041 ) . but increases to values ar ound and above 400 for Prjvf somewhere 
between 0.2 and 0.1 (|Schekochihin et al.ll2005h . There is now also evidence that i? m ,crit 
may actually show a pe ak at Ptm = 0-1 an d might then drop to slightly lower values for 
Pi'M = 0.05 and below (Ilskakovll2007h . This unusual behavior is connected with a change 
of the "roughness" of the velocity field ( Boldyrev fc Cattaneoll2004 ) and the occ urrence 
of a bottleneck effect in the velocity spectrum (jFalkovich 1994; D obler et al.ll2003l ). which 
means that the velocity has maximum roughness for Prjvf ~ 0.1 when the resistive scale 
coincides with the position of the bottleneck. 

The situation is quite different with large-scale dynamos that operate co mpletely inde- 
pendently of the value of Pr^/, as long as ReM is large enough. Already in Brandenburgl 
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Figure 2. Snapshots of B y in the early phase (left ) and satura t ed ph ase (right) 
of a convective dynamo with shear. Adapted from iKapyla et al.l IMPOST ); see also 
http : //www.helsinki . f i/~kapyla/movies .html 



( 20011 ) the critical value of Re m was found to be a round un i ty reg ardle ss of whether 
Ptm = 1 or 0.1. This finding was then extended by Mininni ( 2007 ) and Brandenburg! 
( 2009), who demonstra ted dynamo action d own to PrM = 0.005 and 0.001, respectively, 
or up to Ptm = 1000 ( Brandenburg! [2011 ) . The conclusion is that large-scale dynamo 
action depends solely on the dynamo number, which is given by D = C a for a 2 dynamos 
and by D = C a Cs for afl (or a-shear) dynamos. Here, 



C a = a /rj T ki and C s = S/r/ T kl, 



(2.2) 



where t/t — rj t + n is the sum of turbulent and microscopic magnetic diffusivities, «o is 
a typical value of the a effect, and S is the shear rate (i.e. a typical value of the velocity 
gradient). For C a and Cs, we use standard estimates: 



a 
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where ef = u> • u/k[U^ ms is a measure of the relative kinetic helicity, t 
the turnover time, and 
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With this we find 



where 



|ef»rms/fcl 
"rms/3fcf + T) 



ki 



(2.3) 

(Wrmsfcf) -1 is 

(2.4) 
(2.5) 



i = 1/ (1 + 3Re^) 

is a correction factor that is close to unity for Re a/ 3> 1. Furthermore, we have 

2 / \ 2 

C S 



S/k\ 



Wrms/3fcf + 7/ U ms kf 
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(2.7) 



where we have defined the shear parameter Sh ~ S/u rms k{. Note that, especially in the 
presence of shear, the possibility of dynamo action is strongly connected with the scale 
separation ratio. Indeed, 

(2.8) 



D = -3te f Sh ( ^ 
ki 



depends cubicly on the scale separation ratio. This explains why afl dynamos are often 
much easier to obtain than a 2 dynamos. 
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3. Large-scale dynamos in Cartesian domains 

3.1. Dynamos of aQ type 

Given the alarming signs of earlier investigations by Cattaneo fc Hughes! ( 2006f ) and 
Hughes fc Cattanea (|2008h . it was quite unclear whether the a effect even exists in simu- 
lations with convection. At the time, several po ssible reasons were p ut forward, including 
the absence of stratification; see, for example, Brandenburg! ( 2009 ). In the wake of this 
initial frustration, it was quite surprising wh en large-scale dynam o action was found in 



rotating convection in t he presence of shea r dKapvla, et al" 2008 ); see Figure [2] Similar 



results were obtained bv lHughes fc Proctor (2009J). This controversy was still ongoing at 
the conference "Turbulence and Dynamos" in Stockholm in March 2008^ where Hughe^ 
argued that no convective large-scale dynamos exist, while KapyljQ showed results from 
low Reynolds number convection with shear where large-scale fields were indeed ob- 
tained. In an attempt to clarify the sti ll confli cting results regarding the actual value of 
a in rotating convection, iKapyla, et al.l (|2010aj) pointed out that for a nonuniform mean 
field, the mean current density cannot be neglected. In this case, the turbulent mag- 
netic diffusivity contribu tes and explains the small net electromotive force measured by 
Hughes k, Proctor (2009) by imposing a uniform field. 



3.2. Dynamos of a 2 type 

The simulations mentioned above do not provide conclusive evidence for the existence of 
an a effect in rotating convection, because it is in principle possible that the dynamo could 
be the result of an incoherent att dynamo or a shear-current dynamo. In the absence of 
shear, however, there is no viable alternative to an a 2 dynamo. It is therefore important 
to consider the c o ncept ually simpler case without imposed sh ear, as was also emph asized 
bv lHughes et al. ( 2011 ). who noted that this was not done bv IKapyla et al. ( 2010a ). who 
just considered the cas e of a sinusoidal shear profile. For this re ason, we discu s s in th e 
following the papers of lKapvla etail (|2009al) and, in particular, IKapyla, et al.l I (|2009bl) . 
where large-scale dynamo action was studied in non-shearing convection at sufficiently 
large Coriolis numbers. 

Before trying to simulate an a 2 dynamo for rotating convection, it is instructive to 
obtain guidance from numerically obtained measurements of th e a and turbulent dif- 
fusiv ity tensors. This can be done using the test-field method ( Schrinner et al.l 120051 
20071) , which has been applied t o turbulence in a number of recent papers jBrandenburgl 
2005t i Brandcn burg et al.l 120081 ) . The result is shown in Figure [5] Using this method, 
Kapvla, et al. ( 2009al ) noted that the magnitude of the relevant components of the a 



tensor vary only weakly with Coriolis number, Co = 2f2/u rms fcf, where is angular ve- 
locity, while ?7t diminishes with increasing values of Co. This was a clear indication that 
dynamos of a 2 type might become possible once Co is large enough. We emphasize this 
point, because it is one of the several examples where mean-field theory has proven its 
predictive power. 

Consequently, in a subsequent investigation, IKapyla, et al . (2009b) carried out simula- 
tions for large enough values of Co and did indeed find dynamo action of large-scale type 
when Co il 10. The large-scale field became even more pronounced as the aspect ratio 
was increased. In Figure [3] we present horizontal spectra of magnetic and kinetic ener- 
gies. What is important here is the fact that, even though the magnetic energy is less (by 
factor 5) than the kinetic energy at what we estimate to be fcf (about 5fci), the magnetic 

f http : //agenda . albanova. se/conf erenceDisplay .py?conf Id=325 
J http : //videos .nordita. org/conference/Turbulence2008/hires/Marchl7/Partl . WMV 
% http : //videos .nordita. org/ conf erence/Turbulence2008/hires/March!9/Part5 . WMV 
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Figure 3. Two-dimensional power spectra of 
velocity (upper panel) and magnetic field (lower 
panel) as functions of system size. In the lower 
panel the upper curves show the spectra from 
the saturated state whereas the lower curves 
show the spectra from the ki nematic state 
multip lied by 10 7 . Adapted from iK&pvla, et aD 
(20091). 
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Figure 4. Top panel: rms-values of the to- 
tal magnetic field as functions of time for ver- 
tical field (solid lines) and perfect conductor 
boundary conditions (dashed lines). The two 
lower panels show the sums of the rms-values 
of the Fourier amplitudes of B x and B v for 
k/ki — (middle panel) and k/ki = 1 (bot- 
tom panel). The dotted lines in the two lower 
panels show a saturation predictor according to 
the model of Branden burg (2001). Adapted from 
iKapvla et all (|2009bh . 



energy strongly exceeds the kinetic energy at the scale of the domain. This seems to ex- 
clude alternative explanations whereby the magnetic field at smaller wavenumbers might 
just be a trivial result of diffusion in wavenumber space. Instead, we argue that this is 
strong evidence for the physical reality of an a 2 dynamo driven by rotating convection. 

In agreement with virtually all earlier work on large-scale dynamos o f a 2 type the 
satura tion time of the dynamo is comparable with the resistive time. Indeed, Brandenburg! 
( 200lh found that in the absence of strong magnetic helicity fluxes, the saturation of an 
a 2 dynamo follows a switch-on behavior where, after saturation, the mean field is given 
by 



U 1 

B 2 



1 - e 



-277*4(4-*,,) 



(3-1) 



This is also seen in the present case; see Figure 01 where we overplot the prediction from 
Equation ([53]) . 

It is likely that diff usive magnetic helicity flu xes are present in the convection simula- 
tions discussed above ([Brandenburg et al.ll2009l ) . Those fluxes could in principle give rise 
to faster saturation time s than what is seen in F igure | 3. II This question has bee n ad- 
dressed quantitatively bv lMitra et al.l ( 2010al ) and lHubbard fc Brandenburg! ( 2010t ). who 
note that at the magnetic Reynolds numbers accessible so far, diffusive helicity fluxes 
are still quite weak compared with the resistive processes. Based on their scalings for 
different values of ReM, they estimate that resistive saturation effects would only begin 
to be alleviated for ReM w ell in excess of values around 1000 or even 10 4 ; see also Fig. 10 
of Candelaresi et al. Convection simulations with closed magnetic boundaries do 
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Figure 5. Normalized profiles of a (left panel) and 7/ t (right panel) from kinematic test field 
simulations. The vertical dotted lines at z/d — and z ld = 1 indicate the b ase and top of the 
convectively unstable layer, respectively. Adapted from lKapyla et all (|2009bf ). 



seem t o suffer from catastrophic quenching at least up to ReM ~ 240 ([Kapvla, et al 
2010bl) . Reaching much higher values to study this issue further in the near future is 
not possible. While this is certainly somewhat discouraging news, it does highlight the 
importance of studying detailed scaling properties of large-scale dynamos rather than 
producing solitary examples of dynamos at large resolution, hoping that they represent 
the Sun. 



4. Bi-helical magnetic field 

An important property of a effect dynamos is the fact that they produce bi-helical 
magnetic fields. This means that one expects to see magnetic helicity fluxes with op- 
posite signs at small and large scales. While this is now wel l established theoretically 
(Bla ckman fc Brandenburd2003l:lYousef fc Brandenburg|2003l ). there is still no widespread 
observational evidence for this. Large-scale magnetic helicity can be estimated using syn- 
optic maps of the a zimut hally averaged radial magnetic field of the Sun; see Fig. 1 of 
Brandenburg et al.l (|2003h . Small-scale magnetic helicity fluxes have be en inferred from 



magn etic field measurements in active regions and coronal mass ejections ( Berger fc RuzmaikirJ 
20001) . 



There is also some evidence from me asurements of magnetic h elicity in the solar wind. 
Using the assumption of homogeneity, iMatthaeus et al.f ( 1982h were able to determine 
magnetic helicity spectra for the solar wind. Preliminary analysis of more recent solar 
wind data from the Ulysses spacecraft does indeed suggest that the field in the solar 
wind is also bi-helical. Fu rther details on this are presented in a dedicated publication 
(jBrandenburg et al.|[201lh . 



5. Concluding remarks 

We still do not know exactly how the solar dynamo works. If it is of aVL type, given that 
the a effect is positive in the northern hemisphere, and using the fact that radial shear is 
positive in the bulk of the convection zone, one would expect poleward migration of the 
dynamo wave. This is also what three-dimensional simulations in s pherical shells have 
shown repeatedly over several decades starting with the early work of iGilmanl (|1983l ). and 
now again in the spherical wedge simulations of lKapyla, et al.l (|2010q) . For rapid rotation, 
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however, polarity flips of toroid a l mag netic field can also occur mo re abruptly, as has been 
demonstrated by iBrown et all (|201of) and iGhizaru et~ai1 (|2010|). which is beginni ng to 
be reminiscent of polarity reversals in the geodynamo ( Glatzmaier fc Robertdfl995 ), but 
is different from what we know about the solar dynamo. In this connection it is worth 
recalling yet another recent surprise: oscillatory solutions with equatorward migration are 
also possible in the absence of any differential rotation provide d the dynamo is som ehow 
bounded between highly conducting media at high latitudes (jMitra et al.l l2010bf) . It is 
obviously unclear whether this has any bearing on the solar dynamo, but it reminds us 
of the possibility of surprises. 

Other proposals for how the solar dynamo might work include the flux transport 
dynamo (|Durnevlll995 : Choudhuri et al. 19951 Dikpati fc Charbonneau 19991 ). and the 
possibility of a dynamo shaped by the nega tive radial angular velocity gradient in the 
near-surface shear layer (jBrandenburd 120051 ) . Neither of these two scenarios have been 
seen in three-dimensional turbulence simulations. The former suffers from the difficulty 
of obtaining a sufficiently coherent meridional circulation that does not break up into 
smaller circulation patterns, while the latter may suffer from the difficulty of resolving 
the small-scale turbulence in the near-surface shear layer. A possible step forward might 
therefore be a combined effort utilizing a range of different simulations in spherical and 
Cartesian geometries on the one hand, and improved mean-field theory on the other. 
Clearly, in order to improve mean-field theory it is essential to seek guidance from direct 
simulations, as has already been done in recent years with increasing success. 
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Discussion 

M. PROCTOR: What do you do about the gauge when calculating magnetic helicity, 
and does it make any difference to the answers obtained? 

A. Brandenburg: Yes, the magnetic helicity density is in general gauge-dependent. 
However, if there is sufficient scale separation between mean and fluctuating fields, the 
magnetic helicity density computed from the fluctuating field can be shown to b e gauge- 
invariant ( Subramanian fc Brandenburgll2006 ). Hubbard fc Brandenburg! ( 2010h have re- 



cently confirmed this in a simulation where the magnetic helicity from the mean field 
was strongly gauge dependent, but that from the fluctuating field was not. 

T. Rogers: To calculate magnetic helicity you assumed that the solar wind was 
isotropic, which observations show it is not. How would this affect the results you present? 

A. Brandenburg: Since we have to adopt the Taylor hypothesis, only the two magnetic 
field components perpendicular to the radial direction enter the calculation. The field in 
the plane perpendicular to the radial direction is still fairly isotropic, so I guess our results 
are still meaningful. To clarify the significance of our results further, it might be useful 
to compute magnetic helicity from simulations of anisotropic MHD turbulence with one 
preferred direction. 

D. Hughes: In your simulations that show the generation of large-scale fields on a long 
time, what is the timescale for the generation of the fields? If it is ohmic then it is not 
surprising. 

A. Brandenburg: The initial exponential growth occurs always on a dynamical time 
scale, but full saturation is only obtained on a resistive time scale. We know that this 
problem can only be alleviated by magnetic helicity fluxes, which are quite weak in our 
Cartesian simulations. Nevertheless, our simulations prove the point that the a effect 
works in rotating convection, which was until now quite unclear. 

C. Forest: How do you model the boundary conditions? Open or Closed? 

A. Brandenburg: At the bottom we adopt a perfect conductor boundary condition and 
at the top we assume that the horizontal field vanishes. This pseudo-vacuum condition 
is numerically more robust than a proper vacuum condition. However, it would be more 
realistic to couple the convection simulation to a force-free model, as has been d one for 
forced turbulence simulations in the paper by IWarnecke fc Brandenburg! (|2010l) , which 
is also presented here as a poster. 



